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A search for the lepton flavour violating decays r~ — > l~K% (I = e or y) has been performed 
using a data sample corresponding to an integrated luminosity of 469 fb _1 , collected with the BABAR 
detector at the SLAC PEP-II e + e~ asymmetric energy collider. No statistically significant signal 
has been observed in either channel and the estimated upper limits on branching fractions are 
B{t~ -> e-K° s ) < 3.3 x 1CT S and B(t~ -> yT K° s ) < 4.0 x 1CT 8 at 90% confidence level. 

PACS numbers: 11.30.Fs; 13.35.Dx; 14.60.Fg. 



In the Standard Model (SM), lepton-flavor- violating (LFV) decays of charged leptons are forbidden or highly 
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suppressed even if neutrino mixing is taken into ac- 
count P, @, Any occurrences of LFV decays with 
measurable branching fractions (BFs) would be a clear 
sign of new physics. No signal has been found in exten- 
sive searches for LFV in /i and r decays (e.g. fj, — > 
r — > /X7 d, @, 0]). However, within the bounds set 
by searches, some physics models that extend the SM 
include new sizable LFV processes. For a review, see 
Ref. Q. In this paper a search for r~ — » l~ K® decays is 
presented 



The 



I K~ BF has been estimated in SM ex- 



tensions with heavy singlet Dirac neutrinos 1 1(1 and in 
i?-parity violating supersymmetric models [111 ]. In the 
first case, heavy neutrinos with large mass and large 
mixing with SM leptons are introduced. Because of the 
large number of independent angles and phases in the 
enlarged mixing matrix, the LFV amplitude cannot be 
precisely evaluated. In the large-mass limit of heavy 
neutrinos and keeping only the leading terms, theoret- 
ical upper bound estimations are of the order 10 -16 and 
are thus out of experimental reach. In the second case, 
couplings of SM leptons to new particles are described 
using an i?-parity violating superpotential. With many 
new complex couplings, the phenomenology is immensely 
richer, but at the same time less predictive. While R- 
parity conserving couplings can affect low-energy pro- 
cesses only through loops, i?-parity violating contribu- 
tions can appear as tree-level slepton or squark medi- 
ated processes, competing with SM contributions. So, 
while LFV decays are highly suppressed in the SM, they 
can be significantly enhanced in i?-parity violating su- 
persymmetry. The previous best experimental upper 
limits (ULs) for t~ — > l~K® decay branching ratios 
were measured by the Belle Collaboration using a 281 
fb -1 data sample: B(t~ -> e~Kf) < 5.6 x 10" 8 and 
B(t~ -> H~Kf} < 4.9xl(T 8 at 90% confidence level \r^. 

The measurement described in this paper is performed 
using data collected by the BABAR detector at the PEP- 
II asymmetric energy storage ring. Charged particles are 
detected and their momenta measured by a combination 
of a silicon vertex tracker (SVT) , consisting of 5 layers of 
double-sided detectors, and a 40-layer central drift cham- 
ber (DCH), both operating in a 1.5 T axial magnetic 
field. Charged particle identification (PID) is provided 
by the energy loss in the tracking devices and by the mea- 
sured Cherenkov angle from an internally reflecting ring- 
imaging Cherenkov detector (DIRC) covering the central 
region. Photons are measured, and electrons detected, by 
a CsI(Tl) electromagnetic calorimeter (EMC). The EMC 
is surrounded by an instrumented flux return (IFR) . The 
detector is described in detail elsewhere [l3|, Q ■ 

The analyzed data sample corresponds to an integrated 
luminosity of 469 fb -1 collected from e + e~ collisions, 
425 fb -1 at the T(AS) resonance and 44 fb -1 at center- 
of-mass (CM) energy 10.54 GeV. The total number of 
produced r pairs N TT is (4.31 ± 0.03) x 10 s , calculated 



using the average r cross section of 0.919 ± 0.003 nb esti- 
mated with KK2f [IBj]. The Monte Carlo (MC) simulated 
samples of r leptons are produced using the KK2f gener- 
ator 



in. Decays of 

B mesons are simulated with the EvtGen generator [2(ij ]. 
while e + e~ — > qq events, where q — u,d,s quarks (re- 
ferred to as uds events) or q = c quark, are simulated 
with the JETSET generator [21(. The BaBAR detector 



lil ItJ and Tauola decay library 



is modeled in detail using the GEANT4 simulation pack- 
age 22| . Radiative corrections for signal and background 



processes are simulated using Photos [23]. In the follow- 
ing, the simulated signal and background samples will be 
referred to as signal MC and background MC samples, 
respectively. 

For this analysis, two different stages of selection are 
used. In the first, which we call the loose selection stage, 
we retain enough data to estimate background distribu- 
tion shapes. The second, which we refer to as the tight 
selection, uses criteria that have been chosen to optimize 
the sensitivity. The sensitivity, or expected UL, is de- 
fined as the UL value obtained using the background ex- 
pected from MC: we choose selection criteria that give 
the smallest expected UL. We use loose and tight electron 
and muon PID selectors for the two stages of selection. 
The selectors are based on combinations of measurements 
from the various subdetectors. The average efficiency for 
the loose electron (muon) selector is 98% (92%) for a lab- 
oratory momentum plab > 0.6 (1.4) GeV/c, whereas the 
7r misidentification rate is less than 10% (6%). The aver- 
age identification efficiency for the tight electron (muon) 
selector with a likelihood based algorithm is 93% (80%) 
for the same momentum range, whereas the it misidenti- 
fication rate is less than 0.1% (2%). All selection criteria 
are applied to both channels and quantities are defined 
in the CM system, unless stated otherwise. 

Events are first selected using global event proper- 
ties in order to reject bb, cc and uds background events 
with high multiplicity. All tracks (photons) are required 
to be reconstructed within a fiducial region defined by 
0.410 < 9 < 2.540 (0.410 < 6 < 2.409) radians, where 6 
is the polar angle in the laboratory system with respect 
to the z axis direction [l3[. The overall event charge 
must be zero. Furthermore, the event must include a K® 
candidate with an invariant mass within 25MeV/c 2 of 
the nominal K® mass 2J] , reconstructed from two oppo- 
sitely charged tracks, assuming the pion mass for both. 
The highest momentum track in the CM frame has to 
have a momentum between 1.5 and 4.8GeV/c for both 
modes. For the electron channel events, the total EMC 
energy associated with tracks in the laboratory frame has 
to be less than 9 GeV. The thrust [25} is calculated using 
tracks and calorimeter energy deposits without an asso- 
ciated charged particle track. The thrust magnitude has 
to be between 0.85 (0.88) and 0.98 (0.97) for the elec- 
tron (muon) channel. For each event, two hemispheres 
are defined in the CM frame using the plane perpendicu- 
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lar to the thrust axis. The hemisphere that contains the 
reconstructed r candidate, defined below, is referred to 
as the signal side and the other hemisphere as the tag 
side. Candidate r pair events are required to have three 
reconstructed charged particle tracks on the signal side. 
On the tag side, one track only is required for the muon 
channel, while for the electron channel, events with one 
or three reconstructed tracks are retained. 

The signal r candidates are reconstructed by combin- 
ing one K ° candidate with the third track of the signal 
hemisphere, to which mass is assigned according to the 
considered decay mode. The lepton track is required to 
be identified as an electron or muon by the loose PID 
selector. The signal r candidates are then examined 
in the two dimensional distribution of AE T vs. AM T , 
where AM T is defined as the difference between the in- 
variant mass of the reconstructed r and the world average 
value 24j, and AE T is defined as the difference between 
the energy of the reconstructed r and the expected r en- 
ergy, half the CM total energy. Only r candidates with a 
AM T value within ±0.35 GeV/c 2 and a AE T value within 
±0.4 GeV are retained. The whole decay tree is then fit- 
ted requiring that, within reconstruction uncertainties, 
the K ° decay products form a vertex, the Kg mass is con- 
strained to the nominal value, and the track and the K® 
trajectory form a vertex close to the beam interaction re- 
gion. To improve the energy resolution, a bremsstrahlung 
recovery procedure is applied for the r~ — > e~ K® decay 
mode only: before the fit, the e~ track candidate is com- 
bined with up to three photons with an energy larger 
than 30MeV and contained in a cone around the track 
direction of A6 x Acp = 0.035 x 0.050 rad 2 , where 6 is the 
polar angle and (f> the azimuthal angle in the laboratory 
system. The constrained fit must have a \ 2 probabil- 
ity larger than 1%. If more than one candidate is found 
(which occurs in less than 1% of the events), only that 
with the largest % 2 probability is retained. 

After the above selection is applied, backgrounds re- 
main, mainly from Bhabha events for the electron chan- 
nel and from non-lepton events for the muon channel 
due to the larger pion to muon misidentification. To 
improve the background rejection, further requirements 
are imposed on the K® candidates. For the muon chan- 
nel, the K® laboratory momentum must be greater than 
l.OGeV/c. For the electron channel, in order to remove 
events with a photon conversion faking a if", the in- 
variant mass of the K® daughters, calculated using the 
momentum from the fit and assigning them the electron 
mass, is required to be greater than 0.10 GeV/ c 2 . The 
flight length significance is computed as the three- 
dimensional distance in the laboratory system between 
the r vertex and the K® vertex, divided by its error, and 
we select events with a flight length significance greater 
than 3.0. Finally, the reconstructed mass is required 
to be between 0.482 and 0.514 GeV/c 2 . The last two cri- 
teria are included in the loose selection for the electron 



channel while, for the muon channel, they are applied at 
a later stage in order to maintain sufficient statistics in 
the loose selection sample. The amount of background 
events due to dimuon and Bhabha processes is negligi- 
ble after the loose selection has been applied and most 
of the surviving events come from charm decays, such as 



D~ 



K s n 



and D 



K®£ v. and from combinations 



in the uds events of a true K$ and a fake lepton. 

To avoid bias from adapting selection requirements to 
the data, the tight selection has been optimized in a blind 
way, without looking at the data in the rectangular region 
(blinded box) shown in Fig. [1] corresponding to more 
than ±5 times the resolution for signal events on AE T 
and AM T , respectively. As discussed above, selection 
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FIG. 1: Candidate distributions for signal MC samples r — > 
e~K° (top) and t~ -> yT K% (bottom) in the (AE T , AM T ) 
plane after the loose selection. The rectangle corresponds to 
the blinded box. The z-axis scale is logarithmic. 

criteria have been chosen to optimize the sensitivity on 
the upper limit. Therefore, for the tight selection the 
tighter PID selectors plus the following requirements are 
applied. The event's missing momentum is computed 
by subtracting from the e + e~ momentum all track can- 
didates and all unmatched calorimeter energy deposits. 
To reject events with tracks and photons lost out of the 
acceptance, the missing momentum is required to have a 
transverse component greater than 0.1 (0.2) GeV/c for the 
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electron (muon) channel and the cosine of its polar angle 
in the laboratory system must be smaller than 0.95. In a 
r pair event, when neglecting radiation, the tag-side r has 
the same momentum as the signal-side r but the opposite 
direction. In addition, assuming that the tag r decays to 
a one neutrino (hadronic) mode, the event's missing mo- 
mentum corresponds to the neutrino momentum. These 
two assumptions determine the tag r 4- momentum p TAGl 
and the neutrino 4- momentum p v , respectively, and we 
define the squared invariant mass m^, AG as (p TAG — Pu) 2 ■ 
As shown in Fig. [21 m^ AG peaks at small values for sig- 
nal events and extends to higher values for background 
events. The tail on the right for the signal sample is due 




m 2 AG ((GeV/c 2 ) 2 ) 

FIG. 2: Distributions of m% KG after the loose selection for 
the t~ — > e~ K% channel. The data distribution is shown 
by solid circles with error bars, background MC with a filled 
histogram and signal MC with a dashed line. The signal MC 
distribution is normalized arbitrarily, while the background 
MC is normalized to the data luminosity. The vertical dashed 
line and the arrow indicate the applied requirement. 

to tag t decays to (leptonic) modes with two neutrinos, 
while the tail on the left for the background sample is 
due to events with missing energy from lost photons or 
tracks. The variable rn% AG is required to be smaller than 
2.6 (GeV/c 2 ) 2 for both channels. Shapes for data and 
MC agree within error but a discrepancy is observed in 
the normalization. This does not affect the results be- 
cause the final number of background events is obtained 
using the data sample. The uds background events are 
further reduced by requiring less than six photons on the 
tag side. Signal events have missing momentum due only 
to the undetected neutrino(s) from the tagging r decay. 
Therefore, only for the — » \i~ K® channel, the cosine of 
the angle between the missing momentum and the signal 
r candidate is required to be negative, to further reject 
non-leptonic backgrounds and improve the sensitivity. 

For the final step of analysis, we define another dis- 
criminating variable, Xfull j as the x 2 of the geometrical 
and kinematical fit for the whole decay tree, with addi- 
tional constraints of AM T and AE T equal to 0. Most 
signal events have Xfull values in the range 0-50, and we 



consider this range in the following. In Fig. [3] we show 
the distributions of Xfull f° r data and signal MC inside 
the blinded box after the tight selection. An analytic 
curve describing the background, as detailed in the fol- 
lowing, is also presented. The overall efficiency e in this 
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3: Distributions of Xfull after the tight selection for 
— > e~ K% (top) and t~ — > [i~ K% (bottom) channel. 



The data events are shown by solid circles with error bars. 
The signal MC distributions are shown by dashed lines, while 
the background shapes are shown with filled histograms. The 
signal and background MC distributions are normalized arbi- 
trarily. 

range of Xfull i after the tight selection, and inside the 
blinded box is 9.4% for the t~ -> erK° s mode and 7.0% 
for the — > fi~ K® mode. The total signal efficiency 
is estimated by dividing the number of selected signal 
MC events by the total number of generated r~ — > l~~ K® 
decays and includes the K® — > tt + tt~ BF. 

We estimate the number of background events in the 
signal region using the number of MC background events 
in the range 0-50 of Xfull after the loose selection multi- 
plied by the ratio of numbers of MC background events 
after tight and loose selections in the full range of Xfull- 
We apply a 10% correction to normalize the MC to the 
levels of background seen in data outside the blinded box 
after the tight selection. Total backgrounds of 1.0 ± 0.4 
and 5.3 ± 2.2 events are expected for t~ — > e~K® and 
t~ — * H~K® respectively. Finally the signal region is 
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unblinded and 1 and 2 events are found for electron and 
muon modes, respectively, as already shown in Fig. [3] 

Since no excess above the expected background level 
is found, 90% confidence level limits have been deter- 
mined according to the modified frequentist analysis (or 
method) [26l . \27 \. This method is more powerful 



CL S 

than a simple UL estimation based on numbers of ob- 
served and expected events as it takes into account the 
different distributions of one or more discriminating vari- 
ables between signal and background. The discriminat- 



ing variable used in this analysis is Xfull- The signal 
Xfull distribution is simply provided by the MC sam- 
ple as already shown in Fig. 02 but this cannot be done 
for the background as too few events survive the tight 
selection, but also the loose one. Therefore we obtain 
smooth background shapes by fitting the product of a 
Landau function and a straight line to the MC back- 
ground distributions after the loose selection. Any dis- 
tortions on the shapes that could be introduced by the 
tight selection are negligible compared to the uncertain- 
ties of the shapes themselves. The resulting curves are 
presented in Fig. [3] The adopted test-statistic is the 
likelihood ratio Q = C(S + B)/C{B), where C(B) and 
C(S + B) are, respectively, the likelihood to find the ob- 
served events in the hypothesis of background only and 
of background plus a given amount of signal. The latter, 
and consequently Q, are functions of the hypothesized 
signal BF. The confidence level CL S is defined as the ra- 
tio CL S+B /CL B , where CL S+B and CL B are estimated 
using an ensemble of simulated datasets, generated from 
signal plus background or background only. The genera- 
tion is iterated with a varying hypothetical value of the 
number of signal events, depending on the BF. CL S+B 
and CL B are then the probabilities that the test-statistic 
would be less than the Q exp values observed in data, un- 
der the respective hypothesis. Signal hypotheses corre- 
sponding to CL S < a are rejected at the 1 — a confidence 
level. This method avoids that a negative fluctuation of 
the background is translated into a large improvement of 
the exclusion limit and allows to include uncertainties di- 
rectly on signal and background distributions. The ULs 
on BFs at 90% confidence level are calculated as 



S90 

2eN TT 



(1) 



where sqo is the limit for the signal yield at 90% con- 
fidence level, and e and N TT are already defined above. 
The dominant systematic uncertainties on the signal ef- 
ficiency for the electron (muon) channel come from pos- 
sible data/MC differences in the efficiency of the PID re- 
quirements, 0.4% (5.1%) and of the tracking reconstruc- 
tion, 1.7% (1.6%). Other sources of systematic uncer- 
tainty for the efficiency are: data/MC differences in K® 
reconstruction efficiency (1.0%), the beam energy scale 
and the energy spread (less than 0.2%). The efficiency 
errors from MC statistics are negligible compared with 



the systematics ones. The uncertainty for the total num- 
ber of r pairs comes from the error on the luminosity 
and on the r cross section values (0.7%). We assume 
these uncertainties are uncorrelated and combine them in 
quadrature to give a total signal uncertainty of 2.1% and 
5.5% respectively for the electron and muon channels. 
For each bin of the signal Xfull distribution, we con- 
sider the total uncertainties on the signal yield, and for 
the background distributions the uncertainties on the ex- 
pected background levels. The uncertainties are treated 
as fully correlated between the bins as they are mainly 
due to normalization uncertainties. The analysis results 
are summarized in Fig.|4]presenting CL S for the observed 
events versus the BFs, with the horizontal line defining 
the UL at 90% confidence level. 
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FIG. 4: Observed (full line) and expected (dashed line) CL S 
as a function of the BFs (10 -8 ) for the decays r~ — » e 
and t~ — > (i~ Kg. 



K° s 



From Fig. Gfl the ULs on the BFs at 90% confidence 
level are determined to be: B(t~ — > e~K s r ) < 3.3 x 10~ 8 
and B(t~ -» < 4.0 x 10~ 8 . The CL S obtained 

using the number of expected background MC events, 
instead of data, are shown in the same figure and the 
BF values at 90% confidence level can be regarded as 
the sensitivities: 3.0 x 10~ 8 for the electron channel and 
4.8 x 10~ 8 for the muon one. 

ULs are also determined by exploiting another tech- 
nique that gives a similar but worse sensitivity for the 
UL, so it is used only as cross-check. For this method, se- 
lection criteria on the same quantities were slightly tight- 
ened to reduce the background as much as possible, and 
signal candidates are counted inside the elliptical region 
shown in Fig. [5l The final signal efficiencies with these 
selections are 9.1% for the t~ — > e^K® mode and 6.1% 
for the t~ — » (i~K® mode. The level of background 
in the signal ellipse is estimated by extrapolating the 
event densities found in two sideband regions of AM T , 
as defined in Fig. [5] The AM T background distribution 
is modeled as a linear function plus a Gaussian func- 
tion to account for the peak related to the decay mode 



D~ 



K^tt . This is fitted at the loose selection stage, 
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FIG. 5: Candidate distribution in the (AE T , AM T ) plane 
after all selections for cross-check method (r~ — ► e^ifg on 
the top, r~ — > fi~ K® on the bottom). Data candidates are 
indicated by solid circles. The boxes show the signal MC 
distribution with arbitrary normalization. The blinded box, 
used for both methods, corresponds to the rectangle. The 
gray bands and the ellipse indicate the sidebands used for 
extrapolating the background and the signal region for the 
cross-check measurement. The z-axis scale is linear. 



and B(t- -> (j,-K° s ) < 7.6 x 1(T 8 . 

In conclusion, a search for the lepton flavour violat- 
ing decays r~ — > l~ K® has been performed using a data 
sample of 469fb _1 collected with the BABAR detector at 
the SLAC PEP-II electron-positron storage rings. No 
statistically significant excess of events is observed in ei- 
ther channel and the resulting ULs are B(t~ — > e~K®) < 
3.3 x 1CT 8 and B(t~ -> fi~K°) < 4.0 x lO" 8 at 90% con- 
fidence level. These results are the most restrictive ULs 
on the BFs of these decay modes, and can be used to 
constrain parts of the theoretical phase space in several 
models of physics beyond the Standard Model. 
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where there are sufficient statistics in the sidebands to 
estimate the shape. Then the fitted background distribu- 
tion is normalized according to the number of data events 
in the sidebands after the tight selection. The final esti- 
mated number of background events in the signal region 
is 0.59 ± (0. 19 © 0. 17) and 0.30 ± (0. 17 © 0.05) for the elec- 
tron and muon channels respectively, where the last num- 
ber is the systematic uncertainty accounting for the ob- 
served differences between estimated and real MC sample 
events inside the signal region at the loose selection stage. 
When the signal region is unblinded, we find inside the 
elliptical signal region only one event for each channel. 
Using the signal efficiencies, the estimated residual back- 
grounds, and the number of observed events ULs on the 
BFs at 90% confidence level for this cross-check are cal- 
culated with the POLE program [28| . Uncertainties are 
included assuming that efficiency and background values 
have a Gaussian distribution, and that they are not corre- 
lated. The resulting ULs are B(t~ -> e~K°) < 4.8 xl0~ 8 
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